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Abstract 

We develop a finite temperature field theory formalism in any dimension that has the filling 
fractions as the basic dynamical variables. The formalism efficiently decouples zero temperature 
dynamics from the quantum statistical sums. The zero temperature 'data' is the scattering am- 
plitudes. A saddle point condition leads to an integral equation which is similar in spirit to the 
thermodynamic Bethe ansatz for integrable models, and effectively resums infinite classes of dia- 
grams. We present both relativistic and non-relativistic versions. 
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I. INTRODUCTION 



For the investigation of finite temperature quantum field theory, the standard and well- 
developed formalism is based on the euclidean field theory with time compactified on a 
circle of circumference (3 = 1/T where T is the temperature. In practice, this involves the 
sums over Matsubara frequencies in perturbation theory S]. Though this formalism is 
very useful for some problems, such as the finite temperature dependence of the effective 

n 

potential|3|], for other properties, such as the equation of state involving the pressure and 
density, it would be desirable to have a formalism that more clearly preserves the classical 
picture of a gas of particles at given density and interacting via collisions. 

It helps to realize that at least in principle it is possible to decouple the zero temperature 
dynamics and the quantum statistical sums. The argument is simple: the computation of 
the partition function Z = Tr(e _/3H ) is in principle possible from the complete knowledge of 
the zero temperature eigenstates of the hamiltonian H. In practice this is rather difficult and 
one resorts to perturbative methods such as the Matsubara method, which unfortunately 
entangles the zero temperature dynamics from the quantum statistical mechanics. However 
there does exist a beautiful realization of this kind of decoupling for integrable quantum 
field theories in 2 spacetime dimensions due to Yang and Yang 4|], which is refered to as the 
thermodynamic Bethe ansatz (TBA). (For the relativistic version, see |5|.) In this formalism 
the decoupling is manifest: The free energy is expressed in terms of a pseudo-energy, in fact 
it takes the free field form, and the pseudo-energy is a solution of an integral equation whose 
kernel depends on the zero temperature S-matrix. In fact, the S-matrix and the free-particle 
dispersion relation are the only properties of the theory that are needed as input. The TBA 
is a very powerful tool for tracking RG flows and for computing the conformal central charge 
at the fixed points, and this was one of the main motivations for this work 2s| . 

Dashen, Ma, and Bernstein derived an expression for the partition function in terms of 
the S-matrix in[6]. The derivation is very general, doesn't rely on integrability, and is valid 
in any number of spacetime dimensions 26]. The main result in [6j is rather formal, and a 
considerable amount of additional work is needed to churn it into a useful computational 
iooI. Though some steps toward further developing the formalism were taken in jf| and 
8|, the program was never completed and appears to be eventually abandoned in favor of 
the Matsubara approach. As should be clear in our work, some important aspects of the 
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approach were not well understood in the original works. 

In this paper, we show that, with a new interpretation, the main formula of Dashen et. 
al., can serve as the starting point for a finite temperature formalism based on the physical 
occupation number densities and zero temperature scattering amplitudes. First, by using 
the cluster expansion for the S-matrix in a way not explointed in |6|, we are able to resolve 
certain difficulties previously encountered and render the formalism considerably more ap- 
pealing. Secondly, through a Legendre transformation we are able to formulate the quantum 
statistical mechanics directly in terms of dynamical filling fractions. By "dynamical" we re- 
fer to the property that they are determined from a variational principle, or saddle point 
condition, as in the TBA. 

We w 18 n to also point out the wot, of Lee and Yang on q „ant„ m statists .neehanij 
which preceeded the TBA work of Yang and Yang. That approach is not based on the 
S-matrix but rather on matrix elements of e~^ H itself. Since these matrix elements depend 
on the temperature, that approach does not disentangle the dynamics from the statistical 
sums. Nevertheless, we found some of the ideas, in particular diagramatic description, very 
useful. 

An outline of the paper is the following. In the next section we describe in a completely 
general way how one formulates quantum statistical mechanics of gases in terms of dynamical 
filling fractions. In section III the main result of [6| is reviewed and some potential difficulties 
pointed out with so-called type B terms. In section IV, we show how the extensivity of the 
free energy follows from the cluster decomposition of the S-matrix, however in a delicate way 
that actually provides some constraints on the interpretation. We argue that the formalism 
is only consistent with the cluster decomposition if one passes to euclidean space. This 
also solves the problem with the type B terms. In section V we present a diagrammatic 
description of the formalism. It must be stressed that the resulting diagrams have nothing 
to do with finite temperature Matsubara/Feynman diagrams. In section VI we show that for 
the saddle point construction of section II, one needs only consider the 2-particle irreducible 
diagrams. The saddle point equation is an integral equation that automatically resums 
infinite numbers of diagrams. In section VII we consider contributions that come only 
from 2-body scattering. Here also the integral equation sums infinite numbers of 'foam' 
diagrams. Though we originally had in mind applications to relativistic field theory at high 
temperatures, our construction does not assume the underlying theory is Lorentz invariant. 
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In section VIII we develop the non-relativistic case and solve the integral equation for a 
constant 2-body scattering amplitude. 

In the concluding section we comment on the potential advantages of our formalism for 
certain classes of problems. 

II. FREE ENERGY AS A DYNAMICAL FUNCTIONAL OF FILLING FRAC- 
TIONS 

In the scattering description of quantum statistical mechanics that we will develop, a 
momentum space description is obviously appropriate. A very physical momentum space 
description uses the occupation number densities as the basic dynamical variables. In this 
section we generally describe how this can be done and illustrate it for free particles. 

The free energy density (per volume) T is defined as 

T = --^-\ogZ, Z = Ti e -W~^ (1) 
P * 

where (3 = 1/T, [i are the inverse temperature and chemical potential, V is the d-dimensional 
spacial volume, and H and N are the hamiltonian and particle number operator. Since log Z 
is an extensive quantity, i.e. proportional to the volume, the pressure p of the gas is minus 
the free energy density since p = TdlogZ/dV = —T. For now, let us assume there is one 
species of bosonic (s = 1) or fermionic (s = —1) particle. Given ^(/i), one can compute the 
thermally averaged number density n: 




where k is the d-dimensional momentum vector. (D — d + 1 is the spacetime dimension.) 
The dimensionless quantities / are sometimes called the filling fractions. 

One can express J 7 as a functional of / (and /i) in a meaningful way with a Legendre 
transformation. Define 

G = T{n) + nn (3) 

Treating / and fi as independent variables, then using eq. ([2]) one has that d^G = which 
implies it can be expressed only in terms of / and satisfies 5G/5f = fi. 

Inverting the above construction shows that there exists a functional F{f,fj) 

F(/,A*)=G(/)-/i /(^/( k ) ( 4 ) 
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which satisfies eq. (J2J) and is a stationary point with respect to /: 

The above stationary condition is to be viewed as determining / as a function of fi. The 
physical free energy is then T = F evaluated at the solution / to the above equation. We 
will refer to eq. (jSJ) as the saddle point equation since it is suggestive of a saddle point 
approximation to a functional integral: 



Df e- pvF U> « e- pvy: (6) 

Let us illustrate these definitions for a free theory, in a way that will be useful in the 
sequel. In a free theory, the eigenstates of H are multi-particle Fock space states |ki, k 2 ....). 
Let denote the one-particle energy as a function of momentum k. It is well-known that 
the trace over the multi-particle Fock space gives 

™ = | / log ( : - se ~^~^ (7) 

From the definition eq. (j2J) one finds the filling fractions: 



In order to find the functional F{f,fi) one first computes G from eq. (j3J) and eliminates /i 
to express it in terms of / using eq. (JHJ). One finds 

Fo(/,M) = / 7^((wk-M)/-^[(/ + *)log(l + «/)-/log/]) (9) 

One can then easily verify that SF/Sf = has the solution / = / and plugging this back 
into eq. ([9]) gives the correct result eq. (J71) for JF . 

There is another way to view the above construction which involves the entropy. Write 
eq. © as 

F = £-^S (10) 

where £ is the first (u — /x)/ term in eq. (jUJ), which is the energy density, and S is the 
remaining term in brackets. One can show by a standard counting argument, which involves 
the statistics of the particles, that S represents the entropy density of a gas of particles. 
(See for instance 1(3].) 
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In the sequel, it will be convenient to trade the chemical potential variable \x for the 
variable /o : 

M/./o) =-£/|£ (.ktfl + ./) + /log 

In section VI we will express the corrections to F for an interacting theory in terms of 
scattering amplitudes. 




III. FORMAL EXPRESSION FOR Z IN TERMS OF THE S-MATRIX 

The trace that defines Z is computed with respect to a complete set of orthogonal states 
| a): 

(a'\a) = 5 afa , l = ^2\a)(a\ => Z = ^(a| e^^a) (12) 

a a 

Let us separate H into free (H ) and interacting (Hi) parts: 

H = H + H 1 (13) 



Since the states a in eq. (1121) are not required to be eigenstates of H, let us take the trace 
over eigenstates of H : 

H Q \at) = E a \a) (14) 

In the next section we will specialize to plane-wave scattering states, but for the remainder 
of this section one need not be so specific. 

It was shown how to express the thermal trace in terms of the S-matrix in The 
necessary algebraic tools are familiar from the formal theory of scattering 11 , 3]. There is 
a simple derivation of what will turn out to be the essential term which goes as follows. For 
simplicity we first set the chemical potential fi to zero; it can easily be restored at the end 
by letting aj k — > c<j k — //. Define the resolvent operator 

G < B » = E^lTTTe < 15) 

where E is a real variable and e is small and positive. If H is hermitian, then in the limit 

G(E) - G(E) ] = -2m 5{E - H) (16) 
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Assuming that the spectrum of H has E > 0, one then evidently has 

1 f°° 

Z = / dE e~ pE Tr (G(E) - G\E)) (17) 

Henceforth, we will not always display that operators depend on the variable E; all operators 
depend on E except for H, H ,Hi. 

In order to obtain an expression that is meaningful when the trace is taken over the free 
Fock space, one wants to separate out expressions depending on H as much as possible. 
The resolvent satisfies the equation |27T] 

G = G + G H x G (18) 



where G (E) = (E — H + is) 1 . As in standard discussions of scattering (see e.g. [11]). 
define the operator 

T(E) =H 1 + H 1 G(E) H 1 (19) 

It satisfies 

H X G = TG , GHi = G T (20) 

Using eq. (120]) one now has 

G = Go + GoTGo (21) 

and 

i r°° 

Z = Z - -Im dE e- pE Tr (G TG ) (22) 

JO 

where Zo = Tre~^ H ° is the partition function of the free theory. Using now the cyclicity of 
the trace, OeG = —G^, an integration by parts, and ImG = —n5(E — H ), one finds: 

Z = Zo + Z A + Z B (23) 

where 



Z A = -f3Re dEe~ pE Tr (6{E - H )T(E))^ = -pj^ 6 '^" Re ( T aw 



(24) 



The term Zb involves OeT = —TGqT and is thus quadratic in T. The derivation in 6| 
performs further algebraic manipulations of eq. (1221) and elegantly expresses the final result 
for both Za.b in terms of the S-matrix. The construction is summarized in Appendix A, 
where one finds: 

Z = Zo + — J dE e-? E Tr (S^eS) (25) 



where XOeY = X{8eY) — (8eX)Y, and S(E) is an operator valued function of E related 



to the S-matrix 



281 ] . More specifically, 



S(E) = l-2iri5{E - H ) T(E) 
S-\E) = l + 2m5(E-H ) T\E) (26) 

The on-shell matrix elements of S(E) are the usual S-matrix elements: 

(a'\S(E)\a) = (a'\a) - 2m 5^ T a ,, a , iff E = E a (27) 

where o^) a = 5(E a > — E a ). The condition "iff E = E a n is what is referred to as "on-shell". 

A significant amount of work remains in order to obtain a useful calculational tool from 
eq. (1251) . and this will be the main subject of the next section. For now, note that there are 
two types of terms (in addition to Zq) in eq. (I25I) as in eq. (|23|) . where 

1 f°° 

Z A = -- ^ dE e-? E Tr[d E (5(E - H )(T + Tt))] (28) 



J B 



— vk \ dE e~ pE Tr 



((J(E-Ho)Tt) 0e (5(E-H )T) 



We will refer to these as Type A and Type B terms. Because of the S(E — H ) factors, we 
can now trace over eigenstates \a) of Hq and perform the integral over E. Consider first Za- 
Integrating by parts gives 

Z A = -P z Z e ~ PEa Re(T a . a ) + 5^(0)Re((0\T\0)) (29) 

a 

where |0) is the free-particle vacuum of zero energy. Assuming the vacuum is stable, 
(0|T|0) = and the 5^(0) term can be dropped and eq. (EE]) agrees with (El). (5 (B) (0) 
will be regularized below.) 

The type A and B terms are quite different. The type B terms are actually rather peculiar 
since they may potentially spoil the extensivity of the free energy since they don't obviously 
have the same connectivity properties as the type A terms. We will return to this issue in 
the next section and actually propose that they should be discarded. Z B can be simplified 
using the optical theorem, which follows from S^S* = 1: 

T-T ] = -2ni T f 8{E - H ) T (30) 
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Using this, and the cyclicity of the trace, one finds that the terms involving the derivative 
of the (^-functions vanish. Inserting two complete sets of states one then finds: 

Z B = - i7r e~ pEa 5(E a - E a ,) T*,. a d Ba T a , a (31) 

a,a' 

= 2irJ2 e ~ PEa KE a -E a .) Re{T a ,, a )d E Jm{T a ,. a ) 

a,a' 

Note that the Type B terms vanish if IraT = 0. 



IV. EXTENSIVITY OF THE FREE ENERGY AND THE CLUSTER DECOMPO- 
SITION 

In this section we specialize to a trace over plane wave scattering states and describe 
some new features that arise. In particular, in infinite volume, the free energy is expected 
to be extensive, i.e. proportional to the spacial volume V. So, one first must understand 
how all the various types of terms sum up in a way that can be reorganized as the expo- 
nential of something proportional to the volume. Clearly this has to do with properties 
of connectedness. In the present formalism, this property is essentially a consequence of 
the cluster decomposition of the S-matrix. The volume factors will arise from the following 
regularization of the momentum space delta-function: 

(2n) d 5 ( - d \0) = hm J d d * e tx ^- y) = V (32) 

If the expression ( )23l) clusters in the expected way, then the free energy —TlogZ can be 
identified with the sum of all terms with only one power of V. As we will show, requiring 
that the cluster decomposition for (123]) leads to a free energy with this property actually 
provides some constraints on the interpretation of various terms. 



A. Fock space and S-matrix conventions 

Since we are considering a quantum field theory, the Hilbert space of the free theory is a 
Fock space. Let us now fix our normalizations for the free particle states and their scattering 
amplitudes. The creation-annihilation operators satisfy: 

a k al - s 4,a k = (2n) d <f«>(k - k') (33) 



The Hilbert space is then spanned by the multi-particle states 

JO) (34) 



|kik 2 • • • kjv) = (jj y/2u^j 4 X • • • 4. 



satisfying 

' ™ (35) 



#o|kik 2 - • -kjv) = fy^^kjj |kik 2 ---kjv) 



(The factors are a matter of convention.) One has the non-zero inner products: 

AT 

(ki ... 141k! ... k^) = (2n) Nd sP II W - k *) ( 36 ) 

P i=l 

where the sum is over the N\ permutations V of the order of the k^ and p(V) is the degree 
of the permutation such that p — (p — 1) if V involves an even (odd) number of pairwise 
permutations of particles. The above implies the following resolution of the identity: 

°° 1 f 

1 = ]vT / ^JLL^' ' '^K |ki • ■ ■k Af )(k 1 ■ ■ -kAr] (37) 

N=0 ' J 

where for convenience we have defined the notation: 

dk = / — (38) 

(We have chosen our normalization of states so that the above integration measure over k 
is Lorentz invariant. We emphasize however that we are not assuming the theory to be 
Lorentz invariant; a non-relativistic case is worked out in section VIII.) 



B. Cluster decomposition 

In order to simplify the notation, the free particle states |k 1; k 2 , k n ) will be denoted 
as |12...n) and the S-matrix elements as 

(k^ k' 2 , k^SIki, k 2 , km) = 5 , i/ 2 '... n ';12-m (39) 

For the S-matrix the cluster decomposition is based on the physical requirement that parti- 



cles that are causally separated cannot scatter (see for instance [121]) . The cluster decompo- 
sition also ensures that the free energy only depends on connected S-matrix elements. The 
cluster decomposition may be expressed as follows: 

(a'\S\a)= J! s p S c a , im S c a>2 S c a ,. a3 --- (40) 

partitions 
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where the sum is over partitions of the state \a) into clusters \cxi), \cx2), (The number of 

particles in \cti) and {a[\ is not necessarily the same.) The above formula essentially defines 
what is meant by the connected matrix elements S°. The particles are assumed to be stable 
which implies 

(k'^IkO = (kilkx) = (2n) d 2u kl 6 (d \K - kx) = 6™ (41) 
It will be convenient to express the cluster decomposition in terms of T defined by: 



S = 1 + %T 



(42) 



The cluster decomposition for T is then same as for that of S but without the terms involving 
only delta functions which come from the '1' in S — 1 + iT. (One can show that these 
additional terms involving only delta functions are what give rise to Z .) The connected 
matrix elements T c are characterized as having only a single overall momentum and energy 
conserving delta function and cannot be factorized into functions of only a subset of the 
momenta. Below we will need T c in terms of the conventional scattering amplitudes M.\ 



t>;« = ~^ a T a , a = (2n) d+1 5i E U^(K - k Q ,) M a , a 



(43) 



where E a , k a are the total energy and momentum of the state |a) and 8^} = 8^ E \E a i — E a ). 
For 2^2 and 3^3 particles the cluster decomposition for T then reads: 



Tl'2'i 



12 



J 1'2'; 



12 



7l'2'3';123 — ^1'2'3':123 + ^l'l ^2'3';23 ^ P^rm. 



J 9 



(44) 



The subscript 9 indicates the number of permutations within the brace. If the one-particle 
states were not stable, there would be additional terms. In order to illustrate some important 
additional features, we will also need the 4^4 particle decomposition: 



4';1234 — -^1'2'3'4';1234 ^ -^l'2';12-^3'4';34 P^ rm 



+ 



18 



^i'i ^ 1 2'3'4';234 ^ perm. 



+ 



S V1 S V2 ^3'4';34 ± P erm - 



J 72 



J 16 

(45) 



C. One-particle resummation and emergence of filling fractions 



To compute the required trace, one must set {k' = k}. This leads to a more specialized 
cluster decomposition that is suitable to compute Z. For instance, the 9 terms in braces 
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in eq. (jUJ) separate into 3 + 6 distinct types of terms depending on whether they contain 
«JW(0): 



T 



123;123 



n 



23.123 



+ 



(27i) d 5 w (0)2cu kl f< ± perm. + sS^T* ± perm. (46) 



In order to illustrate the variety of terms that can arise from the cluster decomposition, 
let us compute the Type A terms for low numbers of particles. For simplicity we set the 
chemical potential to zero. For two particles one finds: 



2 part. 



pv 

2 



dki dk.2 e 



-/3(u>i+cj 2 ) 



Re(M 



12:12 J 



(47) 



where we have used eq. ( 1321) . (u>i = c<j ki 
For 3 particles one finds 



3 part. 



^ J dkidk2dk3e-P {uJl+UJ2+W3) Re(M 



123;123J 



J3V 
2 



fs^- / dk l( lk 2 e -/9(«i+«a)( e -/*-i + e -^ 2 )Re(M 



12:12, 



12;12j 



-/3u; k 



(48) 
(49) 
1(50) 



(27r) d / 

These terms have the following interpretation. Eq. f j48~l) is a new contribution of the same 
kind as eq. (14"71) . The term (H9l) just modifies the integration measure dk for the (1471) term. 
The combined measure factors are the first terms in the expansion of the free filling fraction 
fo defined in eq. ([8]). This manner in which the filling fraction emerges was also a feature 



of one approach to finite temperature correlation functions in \l4 . |15| |. The last term is 
proportional to V 2 and is thus expected to be a V 2 term in the expression Z = expl—ftVJ 7 ). 
One can easily verify that it is the correct combinatorial product of a term from the free 
contribution JF and another from the two-particle contribution in eq. ( )47l) . 

A simple combinatorial argument shows that one can sum up all the terms that are a 
product of one T c and some factors to obtain 



T = T ( 



o 




d d K f (k r , 

(2Tr) d 2u; kri 



Re (M 



12..N:12..N, 



+ 



(51) 



N=2 " \n=l 

Restoring the chemical potential fi, the dependence on it is only through fo as given in eq. 



-A 



'"or a scalar field interacting with potential V(<fi) = A0 4 /4!, to lowest order -Mi2 ; i2 = 
131 ] and the 2-particle contribution in eq. fl5T|) gives the same result to order A, as the 
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Matsubara approach^. (In the latter approach this arises as a two-loop finite temperature 
Feynman diagram.) The 3 — > 3 particle contribution to eq. (j51|) was also shown to agree 
with the 3- loop result in jl6|. 



D. The need to continue to euclidean space 



There are additional contributions to T that come from terms in the cluster expansion 
that involve more than one T c factor. They first arise at four particles. The terms in 
7i234 ; i234 that give new contributions to T not already included in eq. (I5TI) are the following: 



1234:1234 



rpc rpc 

1 12:12 1 34:34 



+ 1 



bl 13:12^24:34 



+ i 



12 



rpc rpc 

1 12:34 1 34:12 



+ .... 



(52) 



In order to compute their contribution to Za, we start with the expression flMj) in terms of T. 
We then obtain the cluster decomposition of T from eq. (1521) using T a >. a = —2ir8^) a T a >. a . 
In carrying this out, one sees that one can divide both sides of eq. ( 1521) by 8^ E \0) and 
obtain something well-defined. Finally we express the final result in terms of the scattering 
amplitudes M.. The first term in eq. (15"2"1) in this way leads to: 



(3V 2 



Re 



4 part. 



2m5 {E \0) (^J rfki dk 2 e 



12;12 



+ 



(53) 



Since this is proportional to V 2 it must come from the square of the 2-particle term in the 
expansion of Z — (1 + ...) exp(/3V/2 J Re{M.\2-,i2)) ■ In attempting to match this with eq. 
(|53|) . one notices a very interesting phenomenon: the two terms can only be identified after 
one makes the Wick rotation 

€1 - "«S (54) 
One then regularizes the 5^ E \0) as follows: 



(55) 



The above equation can be understood as following from the fact that in the Matsubara 
approach to finite temperature field theory, one goes to euclidean space and compactifies 
the time on a circle of circumference, i.e. volume, f3. The above equation is then simply the 
euclidean time version of eq. (|3"2"|) . 

Implementing the rules discussed in the last paragraph, the two additional sets of terms 
in the cluster decomposition (1521 lead to two new contributions to the free energy. In the 
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resummation of 1-particle terms, the filling fractions /o again emerge. These lead to: 
s{3 fd^ d d k 2 d% /oCkQ ^g /oCka) Rp(M M , 

T - ~y J (2^p (2^p (2^p ^r^r^r^^ 1 ^ 12 ^ 23 ^ (56) 

- 8 J w? n 2 ^;34^(^ 12; 3 4 ^3 4;12 ) + ... (57) 



where k 4 = ki + k 2 — k 3 and u;4 = oj^. 



E. Z B terms: To B or not to B? 



We now return to issue of the type B terms in eq. (13 ip . It is clear from the previous 
results of this section that the property of the extensivity of the free energy arises in a 
delicate way from the cluster decomposition. Since Z = Z + Za + Zb and the Zg terms 
are quadratic in T they could very easily spoil the extensivity. We have already shown how 
Z + Za exponentiates to something proportional to the volume, but the way this comes 
about is quite subtle. 

We now argue that the need to go to euclidean space cures this potential problem in a 
very satisfying way. In trying to perform the match discussed in the last subsection, one 
actually needs a stronger constraint: 



(3 / Re(M 



12:12, 



Re 



2ti6 (e \0) ( / M 



L12;12 



(58) 



I.e. M i2;i2 must be real. It is well-known that the i mag inary parts of scattering amplitudes 



arise through cuts when one is in Minkowski space 
part vanishes. Therefore, the Zb terms are zero. 



131 ] . In euclidean space the imaginary 



V. DIAGRAMMATIC DESCRIPTION 



The contributions to T have a nice diagrammatic description. Let us represent the 
scattering amplitudes -Mi2.. m; i2..n as a vertex with n incoming and m outgoing lines. We 
construct a diagram with no external lines by linking the lines of vertices. The rules for 
computing a contribution to the free energy density T are then the following: 
(i) Assign a factor of the scattering amplitude and an energy conserving delta function 
27i5^,. a A4 a ';a to each vertex as in Figure 1. 
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(ii) Assign a factor / (k)/2ti;k to each line as in Figure 1. 

(iii) Conserve momentum k at each vertex and integrate over every unconstrained momen- 
tum with J d^k/(2n) d . 

(iv) Identify 2n5^ E \0) = (3. 

(v) Divide by the symmetry factor of the diagram, defined as the number of permutations of 
the internal lines that do not change the topology of the graph, including relative positions. 

(vi) For fermions determine the overall sign of the diagram: it has an overall factor of s if 
it arises from an odd permutation of the outgoing state (1'2'...| in the cluster expansion of 
(1'2'...|T|12...). 

(vii) Divide by —j3. 

1 

= 27r ^12..m;12..n I2..m;12..n 




/o(k)/2ofc 



FIG. 1: Diagrammatic ingredients 

The terms in eq. ( 15T|) are then represented by the diagram in Figure 2. The two terms 
in eqns. ( 1561 ISTj) are represented as diagrams in Figure 3. 

2 

3 




FIG. 2: Diagrammatic representation of terms in eq. (I5ip 

The structure of the "ring" diagrams shown in Figure 4 are especially simple. We compute 
the sum of such diagrams in order to illustrate the above rules and also since we will need the 
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coo 



FIG. 3: Diagrammatic representation of the two terms in eq. (156f) and (|57j) respectively. 

result in the next section. Let the diagram in Figure 4 be denoted as ^Viag- ^he symmetry 
factor of this diagram is 2 for n — 1 and n for n > 2. For fermions, the overall sign of the 
diagram is s n+1 . The sum over n then simply gives a log: 



yo(k) = = y (2^f ■ Mi2 ^ i2(k ' k W ( 59 ) 



FIG. 4: Ring diagrams. 

VI. F AND THE SADDLE POINT EQUATIONS 

We now return to including interactions in the free energy functional F of section II. Let 
us write 

F{f,fo) = F (fJa)+Fi (60) 

where Fo is given in eq. (TTTT) and we define U as the "potential" which incorporates inter- 
actions: 

It is not too difficult to understand that Fi is given by the 2-particle irreducible diagrams 
of the last section. These are defined as diagrams that cannot become disconnected by the 
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cutting of 2 internal lines. For instance, starting from the N = 2 diagram in Figure 2, 
it is clear that the ring diagrams can generated by attaching additional loops. The ring 
diagrams are not 2-particle irreducible and should thus not be included in F i- A more 
detailed argument was given by Lee and Yang [9]. We express this explicitly as follows: F i 
is just ^2-part. irred.(/o) with f replaced by /: 

Fl{f) = -^-pait.irred. (/p ~* f) (62) 

This is shown in Figure 4, where it is implicit that the lines have a factor of / rather than 
fo- 




+ 



FIG. 5: Fx as the sum of 2 particle irreducible diagrams. 



Given F, f is determined by the saddle point equation: 



lo S ( — 7 — ) - lo g 



dU 



(63) 



f J °V fo 

Substituting the solution of (1631) back into F, it can be written in a variety of ways, depending 
on how one utilizes the relation (l63il . One useful way is: 

1 r d d k 



T = - 



(3 J (27r) d 

i r d d k 



[s\og(l + sf) + (l-fd f )U] 
1+sf 



(64) 



(3 J (2vr) c 



s log 



+ (l-fd,)U 



1 + 

It is convenient to define a pseudo-energy e as the following parameterization of /: 

/ s 1 



eP £ - s 

Then the saddle point equation and free energy density take the form: 

1 dU 



(65) 



e = uj — /i 



T = -- 



(3df 
d d k 



pJ-^-A-^g(l-se-^) + (l-fd f )U] 



(66) 
(67) 
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VII. TWO-BODY APPROXIMATION 



If the gas is not too dense, one expects that 2-particle scattering will give the most 
important contribution. This is especially true in the non-relativistic case where scattering 
preserves the number of particles. 

Define the following kernel from the 2-particle scattering amplitude: 

1 



and the convolution: 



Kfk, k') 



(K*/)(k) = 



d d k' 



-Ml2;12(k,k') 



(2tt) c 



K(k,k')/(k') 



The 2-particle contribution to U can then be written as: 

I7(k) = |/(k) (K*/)(k) 
The saddle point equation and free energy then take the forms: 



(68) 



(69) 



(70) 



log 



1 + 3/ 
/ 



log 



1 + s/c 
/o 



1 



P J (2n) d 
In terms of the pseudo-energy: 



e = uj — fi — K * 



dk r slog(l + s/) + ^log 



-PK*f 

1 + sf fp 
f 1 + sf 



Pf. 



-s log(l - se~ Pe ) + ^-{e - u + fi) 



(71) 
(72) 

(73) 
(74) 



(3 J (27r) d _ 

In this two-body approximation, the integral equation (1731 resums all diagrams involving 
two-body scattering. These "foam diagrams" are of the kind shown in Figure 6. 

The foam diagrams contain much more than the ring diagrams. As a check, we now show 
how the ring diagrams are contained in the solution of the integral equation. Let us define 
V as log(yg^) = log(l + y) where 

V = ~7~~r~~~T~ (75) 



1 + s/o 

The integral equation (17T|) can then be expanded in powers of y: 



(1 + 2s/ ) 2 

y 7r f — y + 



PK * (/o + (1 + sf )y) 



(76) 
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FIG. 6: Foam diagrams. 



Now let y = yo + yi + ... and plug this into the above equation. Many kinds of terms are 
generated, but to compare with the ring diagrams we only focus on terms of the type in eq. 
(|59|) . To lowest orders one finds 



where y is defined in eq. ( 1591) . Plugging this lowest order solution into eq. (172|) for the free 
energy one finds 



which agrees with the low order expansion of eq. (J59l) . 

VIII. NON-RELATIVISTIC CASE: HARD-CORE BOSONS 

We have not assumed the underlying theory is Lorentz invariant in the above construction. 
However, in the definition eq. of the states |ki • ■ ■ k„) we included factors of \/2lJu in 
order that dk in eq. ( |36i) is Lorentz invariant. This simplifies the comparision with other 
approaches for relativistic models, as we did in section IVC for the 4 theory, since then M. 
are the conventional scattering amplitudes. 



y = y + 2syl + .... 



(77) 




(78) 
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A. Generalities of the non-relativistic case 

For non-relativistic theories, where uo^ = k 2 /2m, it is more conventional to normalize the 

states as 

|k 1 ...k^> = 4 1 ...4j0> (79) 
The formulas of the previous sections still apply but with the modification: 

r d d k r d d k 



We will keep the definition eq (H3l) of M. in terms of T. The formulas for the two-body 
approximation in section VII then all apply but now with: 

K(k,k')=.M 12 , 12 (k,k') (81) 

The above kernel has dimensions of energy • volume which corresponds to an energy dimen- 
sion oil — d (up to velocity factors). 

B. Hard core boson model 

Consider the two-body potential in position space x: 

V(x,x!) = ^(x-x') (82) 
This leads to the second quantized hamiltonian: 

H = j d d x (^W> f • VV> + |V V W) (83) 
To lowest order T = H x and 

(a'lH^a) = -(2n) d 5 {d) (k a - k a ,) M a > ;a (84) 
Expanding the field in terms of annihilation operators: 

^ (X) = liBf e " x ak (85) 

then to lowest order one finds 

Mi 2; i2 = — 7 = K (86) 
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The coupling constant 7 has units of energy x volume. It can be expressed in terms of a 
physical scattering length a as follows. To first order in perturbation theory the differential 
cross-section in the center of mass is: 

da _ m 2 7 2 -_ 3 

dti ~ 4(27r)*-i 1 J 

where k is the magnitude of k for one of the incoming particles. Since a cross section has 
dimensions of length d_1 , we define a such that the cross section is a d_1 when the wavelength 
of the particle is 2n/a: 

a*' 1 (88) 



da 
dtt 

This leads us to make the definition: 



k^2n/a 



-y a d ~ 2 

(89) 



(2vr) d / 2 m 



C. Lowest order solution 



It is interesting to carry out our analysis for arbitrary spacial dimension d > 0. Integrals 
over momenta can be traded for integrals over u: 

d d k fm\d/i 1 f°° 2V2 

/ duj u {d - 2)/2 (90) 







(2yr) d V27T/ r(d/2) 

For a constant kernel K = —7, and d > 0, the solution to the integral equation eq. f!73|) 
takes the simple form: 

e(k) = c^ k - /i + T6 (91) 
where 5 is independent of k and satisfies the equation 29j 



5 = h lA d / 2 {z^Zs) (92) 

where we have defined the fugacities 

Zfl = e^, z s = e- s (93) 

and a renormalized thermal coupling h and thermal wavelength A^: 

\ d-2 ■, 
a \ . 1 
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The function \A v (z) is the standard polylogarithm, defined as the appropriate analytic con- 
tinuation of 



n=l 



The equation fl92|) is thus seen as non-trivial, transcendental equation that determines 5 as 
a function of /z, T, and the coupling h. Given the solution <5(/i) of this equation, using eq. 
f[T4"|) the density can be expressed as 

1 TS 

H ^ = (2n\ 2 T ) d ^ Ud/2 ^ Z8) = T (96) 

For d > one can integrate by parts and obtain the following expressions for the free energy 

T = ~ {2ir\*) d / 2 ( Li ( d + 2 )/ 2 ^ 5 ) + ^ U d/2(^z 5 )^j (97) 

The case of 2d is interesting since the eq. ( |92l) then becomes algebraic due to Lii(z) = 
-log(l - z). 
problems, in 21]. 



'hese formulas will be applied to Bose-Einstein condensation, and other 



IX. CONCLUDING REMARKS 

In this paper we have mainly focused on developing the formalism in a general way for 
both relativistic and non-relativistic theories. In the two-body approximation, the main 
result is summarized in the two formulas eqs. (173 II 741) and is quite straightforward to im- 
plement: one computes, or measures, the 2-body forward-scattering kernel K at zero tem- 
perature, solves the integral equation ( 1731) . and the free energy is given in terms of this 
solution in eq. (174p . For non-constant kernels, one will probably need to solve the equations 
numerically. 

A formalism that efficiently disentangles zero temperature dynamics and quantum sta- 
tistical sums potentially has many applications, and we discuss some of them in these con- 
cluding remarks. 

Our formalism is especially well-suited for studies of the effects of interactions in Bose- 
Einstein condensation, since there the filling fractions play a central role. We have already 



obtained some results on this problem and will publish them elsewhere 
presented there also give new insights on the Riemann hypothesis. 



211 ] . The results 
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For high-energy particle physics, our formalism has the following potential advantages 
over the Matsubara approach. The necessity to renormalize ultra-violet divergences is es- 
sentially a zero temperature problem, however the Matsubara approach again entangles the 
zero-temperature renormalization with the quantum statistical mechanics and the issue of 
scheme-dependent results can sometimes be a problem. In our formalism, any need for 
renormalization is carried out at zero temperature and the S-matrix is expressed in terms of 
physical quantities at zero temperature, such as particle masses, etc. Our integral equation 
could also shed some light on the infra-red problems that are common in the Matsubara 
approach and also re^e specia, —ationsfl fl fl Q Q Q. We emphasise that 
since the diagrams in this paper have no relation to finite temperature Feynman diagrams, 
our type of resummation is entirely different than the resummations carried out in these 
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other works. Some aspects of the relativistic case will be published in 
may help to st udy the known results on the free energy of strongly coupled supersymmetric 
gauge theories [24]. 



Our formalism may also be well suited to studies of the quark-gluon plasma 23J currently 
being studied at RHIC, since many of the hadronic cross-sections are known. 
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XI. APPENDIX A 



In this appendix we complete the steps leading to the result (1251) obtained in |6( . Starting 
from eq. (|22|) . the next step is to try and use the relation Gq — = —2iri 5(E — H ) so 
that one can integrate over E. One way to do this uses the operators Q: 



Q(E) = GGq 1 = 1 + G Q T 
Q-\E) = GqG- 1 = 1-Goffi 



(98) 



As we show below, Q is the operator that relates free-particle states to the in-states of 
scattering theory. Introduce the notation XdY = X(dY) — (dX)Y. Then, using OeGq = 
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— Gq and H\Q = T one can readily show that 

fi-^fi = -2G%T (99) 

In order to deal with the complex conjugation needed in eq. (122]) . for any operator 
X(E) define X*(E) as simply X(E) but with ie replaced by —ie. This definition gives 
G*{E) = G{E)\ VL*(E) = 1 + G*H U and (ft- 1 )* = 1 - G^. One then has 

(Q^YdEQ* = -2{GlfT ] (100) 

The last step is to define 

S{E) = (Q^ 1 )*Q(E) (101) 
Using eq. (122| [991 . and the cyclicity of the trace, one obtains 

Z = Z + — dE e~ pE Tr (s-^eSJ (102) 

J 

We now demonstrate that the so-called on-shell matrix elements of the S-operator are 
the usual S-matrix. Consider first the operator Q. Inserting a complete set of states la') 
one finds: 

n(E a )\a) = \a) +Y] )-—T a , a \a') (103) 

t—f E a - E a i + ie 

a' 

where 

T a ' a — ( a '\T(E)\a), iff E = E a (104) 

The condition E = E a is the on-shell condition. Eq. H103[) are the Lippmann-Schwinger 
equations: 

\a) in = n(E a )\a) (105) 

which arise from attempting to solve H\a) in = E a \a) in for the in-states (which also have 
energy E a ). 

Turning now to the S-operator, using identities given above, one can easily establish 

S(E) = l-27ri5(E-H ) T(E) 
S-\E) = l + 2ni5{E-H ) T\E) (106) 

The latter implies 

(a'\S(E)\a) = (a'\a) -2m5(E a - E a ,) T a , a , iff E = E a (107) 
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The equation f 1 1 2 j) together with the identification (11061) of the S-matrix is the main result 
obtained in {(|. 
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